Characterizing glucokinase variant mechanisms using a multiplexed abundance assay.

BACKGROUND: Amino acid substitutions can perturb protein activity in multiple ways. Understanding their mechanistic basis may pinpoint how residues contribute to protein function. Here, we characterize the mechanisms underlying variant effects in human glucokinase (GCK) variants, building on our previous comprehensive study on GCK variant activity. RESULTS: Using a yeast growth-based assay, we score the abundance of 95% of GCK missense and nonsense variants. When combining the abundance scores with our previously determined activity scores, we find that 43% of hypoactive variants also decrease cellular protein abundance. The low-abundance variants are enriched in the large domain, while residues in the small domain are tolerant to mutations with respect to abundance. Instead, many variants in the small domain perturb GCK conformational dynamics which are essential for appropriate activity. CONCLUSIONS: In this study, we identify residues important for GCK metabolic stability and conformational dynamics. These residues could be targeted to modulate GCK activity, and thereby affect glucose homeostasis.

Recognition of the HLA-C*07:359 allele in Taiwanese individuals.

One nucleotide substitution in codon 264 of HLA-C*07:02:01:01 results in a novel allele, HLA-C*07:359.

Engineering Enhanced Antimicrobial Properties in α-Conotoxin RgIA through D-Type Amino Acid Substitution and Incorporation of Lysine and Leucine Residues.

Antimicrobial peptides (AMPs), acknowledged as host defense peptides, constitute a category of predominant cationic peptides prevalent in diverse life forms. This study explored the antibacterial activity of α-conotoxin RgIA, and to enhance its stability and efficacy, D-amino acid substitution was employed, resulting in the synthesis of nine RgIA mutant analogs. Results revealed that several modified RgIA mutants displayed inhibitory efficacy against various pathogenic bacteria and fungi, including Candida tropicalis and Escherichia coli. Mechanistic investigations elucidated that these polypeptides achieved antibacterial effects through the disruption of bacterial cell membranes. The study further assessed the designed peptides' hemolytic activity, cytotoxicity, and safety. Mutants with antibacterial activity exhibited lower hemolytic activity and cytotoxicity, with Pep 8 demonstrating favorable safety in mice. RgIA mutants incorporating D-amino acids exhibited notable stability and adaptability, sustaining antibacterial properties across diverse environmental conditions. This research underscores the potential of the peptide to advance innovative oral antibiotics, offering a novel approach to address bacterial infections.

Impact of the H274Y Substitution on N1, N4, N5, and N8 Neuraminidase Enzymatic Properties and Expression in Reverse Genetic Influenza A Viruses.

The H274Y substitution (N2 numbering) in neuraminidase (NA) N1 confers oseltamivir resistance to A(H1N1) influenza viruses. This resistance has been associated with reduced N1 expression using transfected cells, but the effect of this substitution on the enzymatic properties and on the expression of other group-1-NA subtypes is unknown. The aim of the present study was to evaluate the antiviral resistance, enzymatic properties, and expression of wild-type (WT) and H274Y-substituted NA for each group-1-NA. To this end, viruses with WT or H274Y-substituted NA (N1pdm09 or avian N4, N5 or N8) were generated by reverse genetics, and for each reverse-genetic virus, antiviral susceptibility, NA affinity (Km), and maximum velocity (Vm) were measured. The enzymatic properties were coupled with NA quantification on concentrated reverse genetic viruses using mass spectrometry. The H274Y-NA substitution resulted in highly reduced inhibition by oseltamivir and normal inhibition by zanamivir and laninamivir. This resistance was associated with a reduced affinity for MUNANA substrate and a conserved Vm in all viruses. NA quantification was not significantly different between viruses carrying WT or H274Y-N1, N4 or N8, but was lower for viruses carrying H274Y-N5 compared to those carrying a WT-N5. In conclusion, the H274Y-NA substitution of different group-1-NAs systematically reduced their affinity for MUNANA substrate without a significant impact on NA Vm. The impact of the H274Y-NA substitution on viral NA expression was different according to the studied NA.

Influence of Amino Acid Substitutions in Capsid Proteins of Coxsackievirus B5 on Free Chlorine and Thermal Inactivation.

The sensitivity of enteroviruses to disinfectants varies among genetically similar variants and coincides with amino acid changes in capsid proteins, although the effect of individual substitutions remains unknown. Here, we employed reverse genetics to investigate how amino acid substitutions in coxsackievirus B5 (CVB5) capsid proteins affect the virus' sensitivity to free chlorine and heat treatment. Of ten amino acid changes observed in CVB5 variants with free chlorine resistance, none significantly reduced the chlorine sensitivity, indicating a minor role of the capsid composition in chlorine sensitivity of CVB5. Conversely, a subset of these amino acid changes located at the C-terminal region of viral protein 1 led to reduced heat sensitivity. Cryo-electron microscopy revealed that these changes affect the assembly of intermediate viral states (altered and empty particles), suggesting that the mechanism for reduced heat sensitivity could be related to improved molecular packing of CVB5, resulting in greater stability or altered dynamics of virus uncoating during infection.

Structural effects of charge destabilization and amino acid substitutions in amyloid fragments of CsgA.

The most studied functional amyloid is the CsgA, major curli subunit protein, which is produced by numerous strains of Enterobacteriaceae. Although CsgA sequences are highly conserved, they exhibit species diversity, which reflects the specific evolutionary and functional adaptability of the major curli subunit. Herein, we performed bioinformatics analyses to uncover the differences in the amyloidogenic properties of the R4 fragments in Escherichia coli and Salmonella enterica and proposed four mutants for more detailed studies: M1, M2, M3, and M4. The mutated sequences were characterized by various experimental techniques, such as circular dichroism, ATR-FTIR, FT-Raman, thioflavin T, transmission electron microscopy and confocal microscopy. Additionally, molecular dynamics simulations were performed to determine the role of buffer ions in the aggregation process. Our results demonstrated that the aggregation kinetics, fibril morphology, and overall structure of the peptide were significantly affected by the positions of charged amino acids within the repeat sequences of CsgA. Notably, substituting glycine with lysine resulted in the formation of distinctive spherically packed globular aggregates. The differences in morphology observed are attributed to the influence of phosphate ions, which disrupt the local electrostatic interaction network of the polypeptide chains. This study provides knowledge on the preferential formation of amyloid fibrils based on charge states within the polypeptide chain.

Fit for Purpose Approach To Evaluate Detection of Amino Acid Substitutions in Shotgun Proteomics.

Amino acid substitutions (AASs) alter proteins from their genome-expected sequences. Accumulation of substitutions in proteins underlies numerous diseases and antibiotic mechanisms. Accurate global detection of AASs and their frequencies is crucial for understanding these mechanisms. Shotgun proteomics provides an untargeted method for measuring AASs but introduces biases when extrapolating from the genome to identify AASs. To characterize these biases, we created a "ground-truth" approach using the similarities betweenEscherichia coli and Salmonella typhimurium to model the complexity of AAS detection. Shotgun proteomics on mixed lysates generated libraries representing ∼100,000 peptide-spectra and 4161 peptide sequences with a single AAS and defined stoichiometry. Identifying S. typhimurium peptide-spectra with only the E. coli genome resulted in 64.1% correctly identified library peptides. Specific AASs exhibit variable identification efficiencies. There was no inherent bias from the stoichiometry of the substitutions. Short peptides and AASs localized near peptide termini had poor identification efficiency. We identify a new class of "scissor substitutions" that gain or lose protease cleavage sites. Scissor substitutions also had poor identification efficiency. This ground-truth AAS library reveals various sources of bias, which will guide the application of shotgun proteomics to validate AAS hypotheses.

Selective pressure mediated by influenza virus M158-66 epitope-specific CD8+T cells promotes accumulation of extra-epitopic amino acid substitutions associated with viral resistance to these T cells.

Influenza viruses are notorious for their capacity to evade host immunity. Not only can they evade recognition by virus-neutralizing antibodies, there is also evidence that they accumulate mutations in epitopes recognized by virus-specific CD8+T cells. In addition, we have shown previously that human influenza A viruses were less well recognized than avian influenza viruses by CD8+T cells directed to the highly conserved, HLA-A*02:01 restricted M158-66 epitope located in the Matrix 1 (M1) protein. Amino acid differences at residues outside the epitope were responsible for the differential recognition, and it was hypothesized that this reflected immune adaptation of human influenza viruses to selective pressure exerted by M158-66-specific CD8+T cells in the human population. In the present study, we tested this hypothesis and investigated if selective pressure exerted by M158-66 epitope-specific CD8+T cells could drive mutations at the extra-epitopic residues in vitro. To this end, isogenic influenza A viruses with the M1 gene of a human or an avian influenza virus were serially passaged in human lung epithelial A549 cells that transgenically express the HLA-A*02:01 molecule or not, in the presence or absence of M158-66 epitope-specific CD8+T cells. Especially in the virus with the M1 gene of an avian influenza virus, variants emerged with mutations at the extra-epitopic residues associated with reduced recognition by M158-66-specific T cells as detected by Next Generation Sequencing. Although the emergence of these variants was observed in the absence of selective pressure exerted by M158-66 epitope-specific CD8+T cells, their proportion was much larger in the presence of this selective pressure.

Unravelling the link between SARS-CoV-2 mutation frequencies, patient comorbidities, and structural dynamics.

Genomic surveillance is crucial for tracking emergence and spread of novel variants of pathogens, such as SARS-CoV-2, to inform public health interventions and to enforce control measures. However, in some settings especially in low- and middle- income counties, where sequencing platforms are limited, only certain patients get to be selected for sequencing surveillance. Here, we show that patients with multiple comorbidities potentially harbour SARS-CoV-2 with higher mutation rates and thus deserve more attention for genomic surveillance. The relationship between the patient comorbidities, and type of amino acid mutations was assessed. Correlation analysis showed that there was a significant tendency for mutations to occur within the ORF1a region for patients with higher number of comorbidities. Frequency analysis of the amino acid substitution within ORF1a showed that nsp3 P822L of the PLpro protease was one of the highest occurring mutations. Using molecular dynamics, we simulated that the P822L mutation in PLpro represents a system with lower Root Mean Square Deviation (RMSD) fluctuations, and consistent Radius of gyration (Rg), Solvent Accessible Surface Area (SASA) values-indicate a much stabler protein than the wildtype. The outcome of this study will help determine the relationship between the clinical status of a patient and the mutations of the infecting SARS-CoV-2 virus.

Network of epistatic interactions in an enzyme active site revealed by large-scale deep mutational scanning.

Cooperative interactions between amino acids are critical for protein function. A genetic reflection of cooperativity is epistasis, which is when a change in the amino acid at one position changes the sequence requirements at another position. To assess epistasis within an enzyme active site, we utilized CTX-M ß-lactamase as a model system. CTX-M hydrolyzes ß-lactam antibiotics to provide antibiotic resistance, allowing a simple functional selection for rapid sorting of modified enzymes. We created all pairwise mutations across 17 active site positions in the ß-lactamase enzyme and quantitated the function of variants against two ß-lactam antibiotics using next-generation sequencing. Context-dependent sequence requirements were determined by comparing the antibiotic resistance function of double mutations across the CTX-M active site to their predicted function based on the constituent single mutations, revealing both positive epistasis (synergistic interactions) and negative epistasis (antagonistic interactions) between amino acid substitutions. The resulting trends demonstrate that positive epistasis is present throughout the active site, that epistasis between residues is mediated through substrate interactions, and that residues more tolerant to substitutions serve as generic compensators which are responsible for many cases of positive epistasis. Additionally, we show that a key catalytic residue (Glu166) is amenable to compensatory mutations, and we characterize one such double mutant (E166Y/N170G) that acts by an altered catalytic mechanism. These findings shed light on the unique biochemical factors that drive epistasis within an enzyme active site and will inform enzyme engineering efforts by bridging the gap between amino acid sequence and catalytic function.

Structure-guided engineering of biased-agonism in the human niacin receptor via single amino acid substitution.

The Hydroxycarboxylic acid receptor 2 (HCA2), also known as the niacin receptor or GPR109A, is a prototypical GPCR that plays a central role in the inhibition of lipolytic and atherogenic activities. Its activation also results in vasodilation that is linked to the side-effect of flushing associated with dyslipidemia drugs such as niacin. GPR109A continues to be a target for developing potential therapeutics in dyslipidemia with minimized flushing response. Here, we present cryo-EM structures of the GPR109A in complex with dyslipidemia drugs, niacin or acipimox, non-flushing agonists, MK6892 or GSK256073, and recently approved psoriasis drug, monomethyl fumarate (MMF). These structures elucidate the binding mechanism of agonists, molecular basis of receptor activation, and insights into biased signaling elicited by some of the agonists. The structural framework also allows us to engineer receptor mutants that exhibit G-protein signaling bias, and therefore, our study may help in structure-guided drug discovery efforts targeting this receptor.

Enhancing prime editor activity by directed protein evolution in yeast.

Prime editing is a highly versatile genome editing technology that enables the introduction of base substitutions, insertions, and deletions. However, compared to traditional Cas9 nucleases prime editors (PEs) are less active. In this study we use OrthoRep, a yeast-based platform for directed protein evolution, to enhance the editing efficiency of PEs. After several rounds of evolution with increased selection pressure, we identify multiple mutations that have a positive effect on PE activity in yeast cells and in biochemical assays. Combining the two most effective mutations - the A259D amino acid substitution in nCas9 and the K445T substitution in M-MLV RT - results in the variant PE_Y18. Delivery of PE_Y18, encoded on DNA, mRNA or as a ribonucleoprotein complex into mammalian cell lines increases editing rates up to 3.5-fold compared to PEmax. In addition, PE_Y18 supports higher prime editing rates when delivered in vivo into the liver or brain. Our study demonstrates proof-of-concept for the application of OrthoRep to optimize genome editing tools in eukaryotic cells.

Construction of an infectious clone for enterovirus A89 and mutagenesis analysis of viral infection and cell binding.

Enterovirus A89 (EV-A89) is an unconventional strain belonging to the Enterovirus A species. Limited research has been conducted on EV-A89, leaving its biological and pathogenic properties unclear. Developing reverse genetic tools for EV-A89 would help to unravel its infection mechanisms and aid in the development of vaccines and anti-viral drugs. In this study, an infectious clone for EV-A89 was successfully constructed and recombinant enterovirus A89 (rEV-A89) was generated. The rEV-A89 exhibited similar characteristics such as growth curve, plaque morphology, and dsRNA expression with parental strain. Four amino acid substitutions were identified in the EV-A89 capsid, which were found to enhance viral infection. Mechanistic studies revealed that these substitutions increased the virus's cell-binding ability. Establishing reverse genetic tools for EV-A89 will significantly contribute to understanding viral infection and developing anti-viral strategies.IMPORTANCEEnterovirus A species contain many human pathogens and have been classified into conventional cluster and unconventional cluster. Most of the research focuses on various conventional members, while understanding of the life cycle and infection characteristics of unconventional viruses is still very limited. In our study, we constructed the infectious cDNA clone and single-round infectious particles for the unconventional EV-A89, allowing us to investigate the biological properties of recombinant viruses. Moreover, we identified key amino acids residues that facilitate EV-A89 infection and elucidate their roles in enhancing viral binding to host cells. The establishment of the reverse genetics system will greatly facilitate future study on the life cycle of EV-A89 and contribute to the development of prophylactic vaccines and anti-viral drugs.

Hb Hebei [α20 (B1) His→Leu; HBA2:C.62A > T]: A Novel Hemoglobin Variant Found during Measurement of Glycated Hemoglobin.

We report a novel hemoglobin (Hb) variant found in a 34-year-old Chinese male during a routine measurement of glycated hemoglobin. The variant resulted in a P3 peak of 27.5% of the total Hb on high performance liquid chromatography (HPLC) with a glycated hemoglobin mode. However, no abnormal Hb peaks were observed in capillary electrophoresis (CE) with 3.1% Hb A2 and 96.9% Hb A. The amino acid substitution was determined by Sanger sequencing as α20 (B1) His→Leu; the corresponding DNA mutation was identified as CAC > CTC at the first position of codon 20 of the α-chain. This is the first description of the mutation, and we have named it Hb Hebei for the region of origin of the proband.

AAV-mediated hepatic expression of SLC30A10 and the Thr95Ile variant attenuates manganese excess and other phenotypes in Slc30a10-deficient mice.

The manganese (Mn) export protein SLC30A10 is essential for Mn excretion via the liver and intestines. Patients with SLC30A10 deficiency develop Mn excess, dystonia, liver disease, and polycythemia. Recent genome-wide association studies revealed a link between the SLC30A10 variant T95I and markers of liver disease. The in vivo relevance of this variant has yet to be investigated. Using in vitro and in vivo models, we explore the impact of the T95I variant on SLC30A10 function. While SLC30A10 I95 expressed at lower levels than T95 in transfected cell lines, both T95 and I95 variants protected cells similarly from Mn-induced toxicity. Adeno-associated virus 8-mediated expression of T95 or I95 SLC30A10 using the liver-specific thyroxine binding globulin promoter normalized liver Mn levels in mice with hepatocyte Slc30a10 deficiency. Furthermore, Adeno-associated virus-mediated expression of T95 or I95 SLC30A10 normalized red blood cell parameters and body weights and attenuated Mn levels and differential gene expression in livers and brains of mice with whole body Slc30a10 deficiency. While our in vivo data do not indicate that the T95I variant significantly compromises SLC30A10 function, it does reinforce the notion that the liver is a key site of SLC30A10 function. It also supports the idea that restoration of hepatic SLC30A10 expression is sufficient to attenuate phenotypes in SLC30A10 deficiency.

Rheostats, toggles, and neutrals, Oh my! A new framework for understanding how amino acid changes modulate protein function.

Advances in personalized medicine and protein engineering require accurately predicting outcomes of amino acid substitutions. Many algorithms correctly predict that evolutionarily-conserved positions show "toggle" substitution phenotypes, which is defined when a few substitutions at that position retain function. In contrast, predictions often fail for substitutions at the less-studied "rheostat" positions, which are defined when different amino acid substitutions at a position sample at least half of the possible functional range. This review describes efforts to understand the impact and significance of rheostat positions: (1) They have been observed in globular soluble, integral membrane, and intrinsically disordered proteins; within single proteins, their prevalence can be up to 40%. (2) Substitutions at rheostat positions can have biological consequences and ∼10% of substitutions gain function. (3) Although both rheostat and "neutral" (defined when all substitutions exhibit wild-type function) positions are nonconserved, the two classes have different evolutionary signatures. (4) Some rheostat positions have pleiotropic effects on function, simultaneously modulating multiple parameters (e.g., altering both affinity and allosteric coupling). (5) In structural studies, substitutions at rheostat positions appear to cause only local perturbations; the overall conformations appear unchanged. (6) Measured functional changes show promising correlations with predicted changes in protein dynamics; the emergent properties of predicted, dynamically coupled amino acid networks might explain some of the complex functional outcomes observed when substituting rheostat positions. Overall, rheostat positions provide unique opportunities for using single substitutions to tune protein function. Future studies of these positions will yield important insights into the protein sequence/function relationship.

Prevalence of persistent SARS-CoV-2 in a large community surveillance study.

Persistent SARS-CoV-2 infections may act as viral reservoirs that could seed future outbreaks1-5, give rise to highly divergent lineages6-8 and contribute to cases with post-acute COVID-19 sequelae (long COVID)9,10. However, the population prevalence of persistent infections, their viral load kinetics and evolutionary dynamics over the course of infections remain largely unknown. Here, using viral sequence data collected as part of a national infection survey, we identified 381 individuals with SARS-CoV-2 RNA at high titre persisting for at least 30 days, of which 54 had viral RNA persisting at least 60 days. We refer to these as 'persistent infections' as available evidence suggests that they represent ongoing viral replication, although the persistence of non-replicating RNA cannot be ruled out in all. Individuals with persistent infection had more than 50% higher odds of self-reporting long COVID than individuals with non-persistent infection. We estimate that 0.1-0.5% of infections may become persistent with typically rebounding high viral loads and last for at least 60 days. In some individuals, we identified many viral amino acid substitutions, indicating periods of strong positive selection, whereas others had no consensus change in the sequences for prolonged periods, consistent with weak selection. Substitutions included mutations that are lineage defining for SARS-CoV-2 variants, at target sites for monoclonal antibodies and/or are commonly found in immunocompromised people11-14. This work has profound implications for understanding and characterizing SARS-CoV-2 infection, epidemiology and evolution.

Estimating amino acid substitution models from genome datasets: a simulation study on the performance of estimated models.

Estimating parameters of amino acid substitution models is a crucial task in bioinformatics. The maximum likelihood (ML) approach has been proposed to estimate amino acid substitution models from large datasets. The quality of newly estimated models is normally assessed by comparing with the existing models in building ML trees. Two important questions remained are the correlation of the estimated models with the true models and the required size of the training datasets to estimate reliable models. In this article, we performed a simulation study to answer these two questions based on simulated data. We simulated genome datasets with different numbers of genes/alignments based on predefined models (called true models) and predefined trees (called true trees). The simulated datasets were used to estimate amino acid substitution model using the ML estimation methods. Our experiments showed that models estimated by the ML methods from simulated datasets with more than 100 genes have high correlations with the true models. The estimated models performed well in building ML trees in comparison with the true models. The results suggest that amino acid substitution models estimated by the ML methods from large genome datasets are a reliable tool for analyzing amino acid sequences.

Substitution Models of Protein Evolution with Selection on Enzymatic Activity.

Substitution models of evolution are necessary for diverse evolutionary analyses including phylogenetic tree and ancestral sequence reconstructions. At the protein level, empirical substitution models are traditionally used due to their simplicity, but they ignore the variability of substitution patterns among protein sites. Next, in order to improve the realism of the modeling of protein evolution, a series of structurally constrained substitution models were presented, but still they usually ignore constraints on the protein activity. Here, we present a substitution model of protein evolution with selection on both protein structure and enzymatic activity, and that can be applied to phylogenetics. In particular, the model considers the binding affinity of the enzyme-substrate complex as well as structural constraints that include the flexibility of structural flaps, hydrogen bonds, amino acids backbone radius of gyration, and solvent-accessible surface area that are quantified through molecular dynamics simulations. We applied the model to the HIV-1 protease and evaluated it by phylogenetic likelihood in comparison with the best-fitting empirical substitution model and a structurally constrained substitution model that ignores the enzymatic activity. We found that accounting for selection on the protein activity improves the fitting of the modeled functional regions with the real observations, especially in data with high molecular identity, which recommends considering constraints on the protein activity in the development of substitution models of evolution.